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A model Hamiltonian that exhibits asymptotic freedom and a bound state, is used to show on example that 
similarity renormalization group procedure can be tuned to improve convergence of perturbative derivation of 
effective Hamiltonians, through adjustment of the generator of the similarity transformation. The improvement 
is measured by comparing the eigenvalues of perturbatively calculated renormalized Hamiltonians that couple 
only a relatively small number of effective basis states, with the exact bound state energy in the model. The 
improved perturbative calculus leads to a few-percent accuracy in a systematic expansion. 
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I. INTRODUCTION 



The matrix of canonical QCD Hamiltonian is too large for a 
direct diagonalization using computers. One cause of the for- 
bidding size is the dynamical coupling among bare degrees 
of freedom across an infinite range of energy scales. This 
difficulty is manifest in all local quantum field theories of 
physical interest through ultraviolet divergences. One possi- 
ble method to attack the divergence problem in the case of 
asymptotically free theories is to first evaluate an effective 
Hamiltonian and then attempt diagonalization of the corre- 
sponding effective matrix. For some well adjusted effective 
degrees of freedom may interact over only a limited range of 
energy scales, whereby a definition of a finite effective eigen- 
value sub-problem becomes conceivable and its subsequent 
solution using computers is not immediately excluded. To de- 
rive the effective Hamiltonians one can use similarity renor- 
malization group (RG) procedure [1]. This article discusses 
the accuracy that the similarity procedure can achieve in eval- 
uating the effective Hamiltonians in perturbation theory. The 
latter is the only systematic method that exists for handling 
the initial canonical QCD Hamiltonian. 

Having a perturbatively calculated effective Hamiltonian, 
?i(A), that strongly couples effective basis states only if their 
free energies, i.e. eigenvalues of certain Tip (A), do not differ 
by more than a finite width parameter A, one can cut out a 
sub-matrix Wx from the matrix of 7i(A) in the effective basis. 
The energy range of the sub-matrix W\ can be limited to few 
or several widths A. Such W\ contains only a fraction of the 
matrix elements of H{X), which is why W\ is called a win- 
dow. One can diagonalize the window numerically and obtain 
the part of the spectrum of 7i(A) that is matched by the se- 
lected window. For example, the energy range of the window 
matrix W\ can be limited to the region where a bound state 
is formed. The question is then how accurate the perturbative 
RG derivation of 7i(A) can be for the matrix elements in such 
Wx. 

The issue was studied previously [2] using a particularly 
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elegant version of the similarity transformation generator, 
which was taken from Wegner's flow equation for Hamilto- 
nians in solid state physics [3]. In the model studied in Ref. 
[2], the Hamiltonian 7i(A) could be calculated numerically 
with arbitrary accuracy. One could also calculate the same 
Hamiltonian in perturbation theory. Then, one could compare 
results of numerical diagonalization of the resulting windows 
cut out from the exact and from the perturbative results for 
W(A). 

The model study [2] showed that Wegner's generator pro- 
duced quite useful first three terms in an effective coupling 
constant expansion for the window Hamiltonian, W\ = Wq + 
9\ W\i + g\ Wx2. Diagonalization of this Wx produced sur- 
prisingly good accuracy on the order of 10%. The bound state 
eigenvalue of Wx matched the exact one with so small error 
even when A was reduced down to the range that corresponded 
to 1 GeV in QCD. Moreover, only a handful of basis states 
were sufficient in the numerical diagonalization. However, 
gx grew in the asymptotically free model when A was made 
small, which is a generic feature for asymptotic freedom, and 
gx became comparable with 1 for A ~ 1 GeV. At the same 
time, the coefficients of the perturbative expansion for matrix 
elements of Wx beyond the first three terms turned out to form 
an alternating series with growing coefficients and, eventually, 
for A still much larger than 1 GeV, results including terms pro- 
portional to g^, gf, etc. turned out to be useless. Higher-than- 
second order terms led to erratic behavior of Wx and its eigen- 
values. This result required a closer inspection to decide if one 
could improve the perturbative part of the similarity approach 
by using other generators than Wegner's, and whether any se- 
rious investment in a vastly more complex QCD calculations 
should ever be made using this method. The problem was not 
only quantitative, with details depending on the model chosen 
for study, but also qualitative. Namely, one needed to know 
if any significant improvement was possible in principle, be- 
cause the growth of gx for decreasing A in asymptotically free 
theories is in basic conflict with the goal of reducing the width 
A down to the binding scale using perturbation theory. Since 
there are not many exactly soluble models with asymptotic 
freedom and bound states, it was important to find out if any 
improvement could be obtained in the model of Ref. [2]. 

The optimization task may be attempted in similarity RG 
procedure by taking advantage of the fact that the generator of 
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the similarity transformation can be chosen in infinitely many 
different ways. It is plausible that the convergence problems 
emerged in the Wegner case because interactions with large 
energy changes were very quickly eliminated and such rapid 
elimination could produce large feedback in the low-energy 
dynamics, too large for easy reproduction through a perturba- 
tive expansion. Therefore, this work is focused on the ques- 
tion if varying the similarity generator, especially the rate at 
which the off-diagonal terms are removed, can improve con- 
vergence of at least the first five terms in the expansion, i.e. all 
terms up to and including order gf. The number 4 is distin- 
guished by the fact that one needs at least four orders of pertur- 
bation theory to simultaneously account for effective masses 
and running coupling in quark-gluon dynamics in QCD. The 
answer found in the model studied here is positive. Namely, 
it is shown below that one can obtain a systematic expansion 
for W\ into powers of g\, and the sum of first 5 terms in the 
expansion produces at least three times better accuracy than in 
Ref. [2]. When an optimized version of the similarity genera- 
tor is identified, one can include higher order terms than 4th. 
For example, it will be shown here that 6th order improves 
the accuracy to 10 times better than achieved previously us- 
ing Wegner's generator in the model. Sensitivity of the effec- 
tive theory to the finite bare cutoff in numerical calculations 
is also reduced by at least an order of magnitude (this will be 
explained later) and control over numerical instabilities that 
emerge in the RG evolution is also enhanced (this aspect is 
only mentioned). The main result reported here, however, is 
not that a class of similarity generators leads to improvement 
in the model chosen for this study, but that the flexibility avail- 
able in the similarity renormalization group procedure does 
lead to considerable options for improved convergence pat- 
tern. 



Since the similarity approach to QCD has already evolved 
in the effective particle version to third order terms [4] and 
systematic studies of bound states will become possible once 
the fourth order terms are brought under numerical control 
[5], the options for improving convergence need to be further 
explored in more complex models. Also, the development of 
Wegner's flow equations in solid state physics and field the- 
ory [6, 7, 8, 9, 10], has already demonstrated that new options 
for optimization may have much broader range of applicabil- 
ity than the QCD bound state problem could suggest by itself. 
Besides, if the window Hamiltonians of interest could be re- 
liably evaluated using perturbation theory, the standard quan- 
tum mechanical calculus could then be employed to describe 
with them the time evolution of selected states with efficiency 
probably hard to match by any other method. 



This paper is organized as follows. Section II describes the 
model. Section III discusses results concerning perturbation 
theory and optimization of the similarity generator, which is 
followed by concluding remarks in Section IV. 



II. MODEL 

Matrix elements of the model studied here are (see Ref. [2] 
for introduction) 



(2.1) 



where £"„ ~b''\h> 1, and n is an integer, with a convention 
that energy equal 1 corresponds to 1 GeV. The model is cut off 
in infrared and ultraviolet regimes by limiting the subscripts, 
M < m,n < N, M being negative and N positive, both of 
much larger magnitude than 1 . The ultraviolet renormalizabil- 
ity of the model, its asymptotic freedom, and its lack of sen- 
sitivity to the infrared cutoff, were discussed in [2]. The bare 
coupling constant, g — 0.060600631, is adjusted to obtain a 
bound state with eigenvalue - 1 GeV with 8 digits of accuracy 
for h = 2, i\/ = —21 and = 16, in exactly the same way as 
in Ref. [2]. With these choices, the Hamiltonian iJ is a 38 x 
38 matrix with 37 positive eigenvalues, and one negative. The 
largest energy scale reaches 65 TeV, and the smallest one is 
0.5 KeV. This energy range saturates the needs of contempo- 
rary theories where asymptotic freedom and bound states are 
of interest. 

The differential similarity RG procedure is based on a class 
of equations that can be written as 



dU 



Fx{n},n 



(2.2) 



with the initial condition set at A = oo, forcing the RG trajec- 
tory to start from the initial Hamiltonian, Ti,{oo) = H, of the 
bare theory with counterterms. The generator of the similarity 
transformation can be written as (!?,„ = Timm) 



fm n{^m ) n 



(2.3) 



where different choices of /,„„ lead to different results for 
the trajectory of effective Hamiltonians. /,„„ = ds/dX, 
where s = 1/A^ and (j>mn = 1, gives Wegner's equation, 
which has a considerable record of applications in solid state 
physics [3, 6, 7, 8, 9], independently of the similarity renor- 
malization group studies of asymptotically free theories. Note 
also, that multiplication of 0m„ by a constant, say a, is equiv- 
alent to change of variables from A to A' = A/ y^, and corre- 
sponds to mere shift on logarithmic scale. 

Behavior of matrix elements 7imn(A) as functions of A has 
already been extensively discussed in the case of Wegner's 
generator in Ref. [2]. For <f)mn that is independent of A, it can 
be qualitatively represented as 
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where gx is the effective coupling constant. The Gaussian- 
like fall-off function is not exact, but it is not far from the 
actual form factor that makes the interaction Hamiltonian ma- 
trix narrow on the energy scale, i.e., of width A. It is clear that 
in the case of Eq. (2.4) the effective coupling constant could 
be evaluated from the following formula, 
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The same definition is applied in the case of exact solution. 
The resulting effective coupling gx can be computed numeri- 
cally in all cases considered here and its dependence on A is 
known, see Fig. 1 . For A 1, one has g\ ^ 1/ log A, and the 





FIG. 2: Bound state energies obtained by diagonalization of exact 
window matrix, Wx, with a) M = -8 N = 2, b) i\/ = -8 iV = 
1, and c) A/ = —8 A'' = 0, as functions of A in case of Wegner's 
generator. 



FIG. 1 : The exact running couphng constant gx from Eq. 2.5 plotted 
as function of the width A. The curve labeled "Wegner" results from 
Wegner's equation and the curve labeled "c = I" is obtained using 
generator defined by Eq. (2.3) with (l>mn given in Eq. 3.4, where 
c= 1. 

growth of g\ toward smaller A continues in case of Wegner's 
generator until the coupling constant slightly exceeds 1 , which 
happens when the effective Hamiltonian width becomes com- 
parable to the momentum-space width of the bound eigenstate 
wave function. For comparison, the optimized generator result 
is labeled c = 1. It will be discussed later 

Near A = 1 GeV the bound state eigenvalue begins to build 
up on the diagonal of H, with participation of a limited num- 
ber of neighboring matrix elements. The limitation comes 
from the size of A and the actual wave function width. Thus, 
when one extracts a sub-matrix from 7i(A ^ 1), which forms 
a window, W\, that embraces all matrix elements that count in 
the bound state dynamics, the diagonalization of that window 
alone can produce the right value of the bound state energy. 
This is shown in Figure 2 in case of Wegner's equation for 
three different choices of the window matrix, i.e. for three 
different choices of AI and N that limit Wx in the same fash- 
ion as M and N limit H in Eq. (2.1). 

Figure 2 demonstrates that an exactly calculated small win- 
dow may provide correct results when A is appropriately 
small, i.e. below about 4 GeV. The only problem is that the 
calculation of a window Hamiltonian with finite A in realis- 
tic theories is not feasible beyond perturbation theory. Thus, 
one is led to the question how well the same small windows 
can be calculated perturbatively. The question is of key im- 
portance to QCD because the region of A < 4 GeV is where 
a significant raise of the effective coupling constant g\ is ex- 
pected. But the question may also be relevant in other areas 
of physics than strong interactions. 



III. PERTURBATION THEORY 

This Section describes difficulties that similarity approach 
encounters with convergence of perturbation theory, treating 
the case with similarity transformation generator taken from 
Wegner's flow equation as a benchmark, and shows how one 
can optimize the convergence by changing the generator The 
quality of the perturbative expansions is established in the 
model by comparing the negative eigenvalue of perturbatively 
calculated windows W\ with the exact value of — 1 GeV. 

Perturbation theory is constructed in the following way. 
First one writes a series for 7Y(A) in the form 



7^(A) = ^g"H(")(A) , 



(3.1) 



where g is the bare coupling constant from H of Eq. (2.1). 
For example, in second order expansion the solution of Eq. 
(2.2) with <f>mn independent of A, reads 

H{X) = E„iS„in 

- g^/EmE„ exp [ - 4>mniE^ - Enf/X^] 



g^^/EmE„e^p [ - 0™„(-B™ - E^f/X^] 

\ik{Em — Ek) ~ 4>kn{Ek — En) 



L k=M 



1 — cxp I 



,-A 



mnk ' 



A2 



{E„i — Enf 



corrections , 



where 
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(3.2) 
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— 4>m7i{E,n ~ EnY , (3.3) 

and corrections are of higher order in g than second. The latter 
terms are also calculable analytically in orders 3rd, 4th, and 
higher, but they are too complex to display here and hardly 
useful in the analytic form for numerical calculations in the 
present study, because they involve combinations of mutually 
canceling terms that contain ratios of functions that approach 
zero in numerator and denominator and require careful eval- 
uation. Instead of such analytic expressions, which are dis- 
cussed here only to demonstrate what the required expansion 
consists in, it is more practical in the model to apply a numer- 
ical algorithm and generate matrices 7i(")(A) order by order 
through Runge-Kutta integration of {N - M + 1){N - M) /2 
coupled non-linear differential equations. 

Once the expansion of Eq. (3.1) is calculated to some or- 
der, one uses Eq. (2.5) to express gx as a series in powers 
of g to the same order and one inverts that series to obtain g 
in the form of a new series in powers of gx up to that same 
order The latter series is then inserted into Eq. (3.1), which 
provides the desired effective coupling constant expansion for 
the whole matrix 7Y(A), i.e. the expansion in terms of pow- 
ers of gx- Note that this procedure can be carried out with- 
out detailed analytic knowledge of how the matrices Ti'") (A) 
depend on A, although such knowledge and selection of rele- 
vant, marginal, and irrelevant matrices may be of great help in 
calculations with more than one coupling constant if one can 
reliably identify correlations among the couplings. In fact, the 
accuracy achieved here in plain perturbation theory can most 
probably be further enhanced using techniques similar to the 
ones employed in Ref. [10]. 

The accuracy of the perturbative procedure for evaluating 
effective Hamiltonians in case of Wegner's equation, i.e. with 
0„j„ = 1, is shown in Fig. 3, in terms of the negative eigenval- 
ues of a selected window Wx ■ While the window matrix ele- 
ments are evaluated in perturbation theory, the eigenvalues of 
the window are obtained from non-perturbative diagonaliza- 
tion. Note that the vertical axis has logarithmic scale. While 
the second order result produces about 10% accuracy after di- 
agonalization, the next orders in the expansion are wrong by 
huge factors and completely unacceptable, including the fact 
that they render spurious negative eigenvalues. Since it is not 
known how to select the right negative eigenvalue of the win- 
dows obtained in third and fourth order calculation, all neg- 
ative eigenvalues of the windows are displayed, labeled with 
the number that indicates the order of the expansion for Wx- 

The window Wx used in Fig. 3 was selected on the basis of 
results analogous to those shown in Fig. 2. An arbitrarily set 
requirement was adopted, that a good window should produce 
the negative eigenvalue with accuracy around 0.5 percent for 
4 > A > 1 when calculated exactly. The 0.5% may be con- 
sidered a large error but it is introduced taking into account 
that low order perturbation theory for Wx is expected to lead 
to much larger errors. The bound state energies calculated by 
diagonalization of exact windows with different values of N 
and one value of M = —8 are shown in Figure 2. It is clear 
that one should choose = 2 to make sure that the correct 
result is contained in the window of choice when A < 4. For 



A'^ kept equal 2 at A = 2, one obtains the eigenvalue —9.9599 
with M = -8, -9.9199 with M = -7, and -9.8398 with 
M = —6. Therefore, M — —8 and N — 2 were selected in 
Fig. 3. 
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FIG. 3: Negative eigenvalues obtained from diagonalization of the 
window Wx with M — —8 and N = 2, obtained with Wegner's 
generator in first four orders of perturbation theory. The orders 3 and 
4 produce spurious negative eigenvalues (see the text for details) with 
no sign of convergence. 

The summary of results obtained for Wx from the pertur- 
bative expansion of Wegner's equation is that the expansion 
produces accuracy on the order of 10% in second order but 
completely fails in higher orders. Therefore, it is important to 
verify what happens when one considers /„i„ with (pmn ^ 1 
in Eq. (2.3). Two basic options can be considered. One is 
to make > 1, which accelerates changes in the matrix 
elements 7irnn(A) as functions of A, and the other is to make 
<l>mn < 1, which slows down the changes. 

Numerical calculations show that the acceleration with a 
factor of the type (/),„„ = 1 + c |to — ri| with c ranging be- 
tween 0.01 and 100 does not change the erratic behavior of 
the window to any better. In contrast, slowing down the evo- 
lution by the factor of the type 



produces significant improvements. The type of results one 
obtains from the slowed down similarity transformation is il- 
lustrated by an example in Fig. 4, using Eq. (3.4) with c = 1. 
In this case, one obtains only one negative eigenvalue in the 
range [—2, 0] (no spurious negative eigenvalues appear in that 
range). Moreover, the accuracy of the eigenvalue calculated 
by diagonalization of the fourth order window reaches 3%, 
which is unattainable using Wegner's generator. 

The particular choice of c = 1 follows from the observation 
that 4th order results are less accurate than the 3rd order for 
c < 0.5, while for c > 2 no further improvement is obtained 
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FIG. 4: Bound state energy obtained by diagonalization of the same 
window Hamiltonian as used in Fig. 3, but with the Hamiltonian 
matrix elements calculated in first six successive orders of perturba- 
tion theory using Eq. (3.4) for the similarity generator, instead of 
Wegner's. 



in the accuracy of the 4th order terms. Figure 4 shows also 
results of including orders 5th and 6th for c — 1. Inclusion 
of 6th order contribution achieves the accuracy of 1 % for the 
bound state eigenvalue, which is 10 times better than with the 
Wegner generator. Since the contribution of 6th order is larger 
than 5th, one may consider larger c, around 2, to obtain a more 
systematic inclusion of successive orders as required. This 
issue will not be further discussed here, since already the 4th 
order of rigorous perturbation theory (two loops) for effective 
Hamiltonians in QCD is currently not known. 

Finally, one has to address the issue of renormalizabiUty 
since for c ^ oo the similarity transformation generator ap- 
proaches zero within the bounds set by the ultraviolet cutoff 
N and no RG evolution occurs. The question is then how 
strongly the effective window Hamiltonians W\ for c ~ 1 
and A < 4 GeV, depend on N. This can be studied by in- 
spection of examples and the result is that the generator with 
c = 1 produces considerably more efficient approach to renor- 
malized effective theory than Wegner's. The measure of how 
quickly the renormalized effective theory is achieved is found 
by evaluating 64 matrix elements of the exact window Wx 
with M = -4 and iV 3 for TV = 16, and for N = 20. 
These matrix elements are most important in the formation of 
bound-state in the effective dynamics. In both cases the bare 
coupling constant is adjusted to produce the same bound state 
eigenvalue of -1 (one obtains g = 0.04878048667 in the case 

~ 20). Then, one evaluates 



1 



{N - M)2 



N 

E 



Wxmn{N = W) 



WxmniN = 20) 



Figs. 3 and 4, the chosen window is most suitable for diag- 
onalization. One cannot rely on a comprison at one value of 
A because, as it is visible in Fig. 1, the formation of bound 
states occurs in Wegner's case (c = 0) and for c 7^ at differ- 
ent values of A, being shifted toward smaller values of A for 
c = 1. It turns out that the resulting measure r(A) for the Weg- 
ner generator varies between 10~^ and 10~^, with numerical 
accuracy being of significance, while for c = 1 one obtains 
a numerically stable result, approximately equal 10~® in the 
whole range between 4 and 1 GeV. 

The utility of logarithmic (i.e. depending on exponents m 
and n of the base 6) slowing down of RG evolution of the off- 
diagonal matrix elements originates from two sources. One 
is the smaller size of the effective coupling constant in the 
region where the bound state is formed, as seen in Fig. 1. For 
example, at A = 2 GeV, in Wegner's case g\ ~ 0.72373722, 
and the slowed down case with c = 1 gives g\ = 0.28518167, 
which is about 2.5 times smaller and produces about 40 times 
smaller result for gf. The other source is that coefficients of 
the perturbative series in powers of gx for matrix elements of 
7i(A) derived with slower generator are not growing with the 
order as they do in Wegner's case. This is illustrated in Tab. 
I on two generic examples. Namely, the matrix elements of 
coefficient matrices a; with i — 0, 1, 2, 3, and 4, in 



niX) =ao + aigx + a2g\ + a^^gl + a^gl + 



(3.6) 



are tabulated for {m,n) equal (—1, —1) and (2, 1) at A = 2 
GeV, in Wegner's case and for c = 1, as indicated. It is visible 
that the ratios a^/a^ are about an order of magnitude smaller 
in size for c = 1 than in Wegner's case. 



A = 2 


F{H} 


ao 


ffli 


02 


as 


a4 




Wegner 


0.5 


-0.498 


0.044 


0.253 


-0.902 


c = 1 


0.5 


-0.499 


-0.689 


-1.417 


-0.462 


7^2,1 


Wegner 





-1.050 


0.270 


6.116 


-24.788 


c = 1 





-1.724 


-2.685 


-2.741 


0.960 



(3.5) 

in the range of A between 4 and 1 GeV where, as indicated by 



TABLE I: Numerical values of perturbative coefficients of successive 
powers of gx in Eq. (3.6) at A = 2 GeV for selected matrix elements 
in case of Wegner's generator and for c = 1. See the text for details. 



IV. CONCLUSION 

Similarity renormalization group approach to solving 
asymptotically free theories could employ various generators 
of the similarity transformation to calculate effective Hamil- 
tonians. The calculation can be carried out in perturbation 
theory and then the Hamiltonians can be diagonalized using 
computers with a reasonable chance for obtaining accurate 
answers, if the perturbative procedure for evaluating effec- 
tive Hamiltonians is sufficiently precise. It was shown here, 
using Wegner's generator as a benchmark, that in order to ob- 
tain stable 3% accuracy in the bound state spectrum one needs 
to slow down the rate of elimination of the off-diagonal ma- 
trix elements. In the model study, the rate of change of bare 
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Hamiltonian matrix elements with largest energy changes was 
slowed down 38 times. The stabilizing effect of slowing down 
the RG flow is evident from comparison of Figs. 3 and 4. Per- 
turbation theory including sixth order terms with c = 1 in Eq. 
(3.4) achieved 1% accuracy in evaluating effective dynamics, 
as measured by the bound state eigenvalue. 

However, there is a stiff price to pay. Namely, the second 
order perturbative expansions with slower generators produce 
only about 30% accuracy, in comparison to 10% in Wegner's 
case, and to take advantage of the improved convergence one 
has to go to higher orders, which requires considerably more 
effort. On the other hand, orders higher than second are not 
attainable using Wegner's generator, and the outstanding ac- 
curacy of about 10%, achievable already in the second or- 



der, cannot be improved by including higher orders in a plain 
expansion in one running coupling constant. Moreover, the 
slower generator with c = 1 leads also to less sensitivity of 
the numerically evaluated effective theory to the bare cutoff, 
when the latter is kept large but finite. 

Detailed studies will disclose how much work is actu- 
ally required to achieve the few-percent accuracy with realis- 
tic Hamiltonians. However, the model study described here 
shows akeady that optimization of the similarity generator 
may in principle lead to a desired convergence pattern in the 
case of bound states in asymptotically free theories. Perhaps, 
the optimization found here may turn out to be useful also in 
other applications of similarity than QCD. 
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